Ideal strength of random alloys from first-principles theory 
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The all-electron exact muffin-tin orbitals method in combination with the coherent-potential ap- 
proximation has been employed to investigate the ideal tensile strengths of elemental V, Mo solids 
and V- and Mo-based random solid solutions. The present ideal tensile strengths, calculated as- 
suming isotropic Poisson contraction, are 16.1, 26.7 and 37.6 GPa for bcc V in the [001], [111] and 
[110] directions, respectively, and 26.7 GPa for bcc Mo in the [001] direction, which are all in good 
agreement with the available theoretical data. When a few percent Tc is introduced in Mo, it is 
found that the ideal strength decreases in the [001] direction. For the V-based alloys, Cr increases 
and Ti decreases the ideal tensile strength in all principal directions. Adding the same concentration 
of Cr and Ti to V leads to ternary alloys with similar ideal strength values as that of pure V. The 
alloying effects on the ideal strength is explained using the electronic band structure. 

PACS numbers: 62.20.-x,71.15.Nc,71.20.Be,71.23.-k 



T3 

C 

O 
o 



> 

CN 
^1" 
O 

CO 

CN 
O 
CO 



I. INTRODUCTION 

High strength and good ductility define the two most 
important mechanical properties of metallic structure 
materials.— In materials, strength is usually controlled by 
the occurrence of grain boundaries, cracks, dislocations 
and other micro-structural defects. If such defects were 
not present, the strength would be limited by the stress 
at which the lattice itself becomes unstable with respect 
to a homogeneous strain. This stress referred to as the 
ideal strength, is the possible maximum strength of an 
ideal single crystal. The ideal strength is an inherent 
property of a material which can offer insight into the 
correlation between the intrinsic chemical bonding and 
the crystal symmetry and has been accepted as an essen- 
tial mechanical parameter of single crystal materials j 2 - 

The experimental data on ideal strength are rather lim- 
ited, since it is difficult to measure the ideal strength us- 
ing experimental tools. There are only few experimental 
data that were obtained from tensile tests for whiskers 3 
and from nanoindentation experiments^ Recently, calcu- 
lations of ideal strength became of great interest because 
they represent an upper bound to the strength of any 
real crystal on the attainable stress. The ideal strength 
may be calculated under special loading conditions: ideal 
strength in tension, e.g., for high symmetry directions of 
cubic crystals under uniaxial load, or ideal strength in 
shear, e.g., for common slip systems < 111 > {112} and 
< 111 > {110} in body-centered cubic (bcc) crystals^ 

The ideal strength is one of few mechanical properties 



that can be calculated from first principles. The ideal 
strength of refractory metals (such as Mo, Nb, V, and 
W), late transition metals (such as Cu, Pt and Au), el- 
ements crystallizing in the diamond structure (Si, Gc, 
and C) as well as some ordered alloys (TiAl and N13AI) 
have been extensively investigated^— For example, Luo 
et a/ji£ and Nagasako et al& focused on the ideal tensile 
strength (ITS) of the bcc metals Mo, Nb, and V in the 
< 001 > directions. Accordingly, in all of them the de- 
formation starts along the Bain path, but branches away 
onto an orthorhombic path before the face-centered cubic 
(fee) point is reached. In Mo, this orthorhombic distor- 
tion does not influence the ITS; however, in Nb and V, 
the branching occurs before the maximum stress on the 
tetragonal deformation curve is reached and hence causes 
a significant decrease in the ITS. The stress-strain rela- 
tions and the corresponding theoretical tensile strengths 
of ordered 7- TiAl (Lin structure) alloy exhibit strong 
anisotropy in different crystalline directions originating 
from the structural anisotropy of 7-TiAlJ^ In spite of all 
these theoretical efforts, the first-principles description 
of the ITS in random solid solutions is rather limited. 
This is because the material strength is sensitive to their 
structure. The only attempt^ employed the virtual crys- 
tal approximation (VCA) to study the ITS of the binary 
Ti-V for concentrations of Ti equal to and higher than 
30 at.%. The results showed that the more the ITS is 
decreased the more Ti is present in the alloy. 

Vanadium rich ternary V-Cr-Ti alloys are important 
candidate structure materials for the first-wall/blanket of 
future fusion reactors i 20 i 21 These V-based alloys exhibit 
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excellent mechanical properties, decent thermal creep be- 
havior, high thermal conductivity, good resistance to 
irradiation-induced swelling and damage, and long op- 
erating lifetime in the fusion environment i 2 ^— Consid- 
erable efforts have been made to find optimal V-based 
alloy compositions that can endure the extreme environ- 
ment of fusion-reactors. The available experimental data 
indicates that reasonable properties can be achieved by 
introducing a few percent Ti and Cr into the V matrix^ 
Consequently, the ternary V-Cr-Ti system has attracted 
broad interest, and in particular the compositions with 
0—15 at.% Cr and — 20 at.% Ti have been intensively 
investigated To the best of our knowledge, most ex- 
perimental efforts on V-Cr-Ti alloys have been devoted 
to the ductile-brittle transition temperature before and 
after irradiation, swelling properties, and impact tough- 
ness as a function of Cr and Ti contents. 

As a promising structure material for fusion reactors, 
V-based alloys must not only withstand radiation dam- 
age but should also keep intrinsic mechanical proper- 
ties and structural strength. Thus it is necessary to 
study fundamental mechanical properties, such as elas- 
tic properties, which were systematically investigated 
for Vi-x-yCr^Tij, random alloys (0 < x < 0.1 and 
0<y<0.1)in our former worksi 32 ' 33 The second- 
order elastic constants describe the mechanical proper- 
ties of materials in the small deformation region, where 
the stress-strain relations are linear. The ideal strength 
describes the mechanical properties of the material be- 
yond the elastic region, which is important during the 
alloy design. 

In this work, we use the all-electron exact muffin- 
tin orbitals method (EMTO) in combination with the 
coherent-potential approximation (CPA) to investigate 
the ITS of bcc elemental V and Mo solids as well as 
Mo . 9 Tc .i and Vi^^Cr^T^ (0 < x + y < 0.1) ran- 
dom alloys as a function of concentrations. For sake of 
simplicity, here we consider only the ideal tensile strength 
corresponding to isotropic Poisson contraction. However, 
the present approach can easily be generalized to account 
for shear instabilities as reported in the case of several 
elemental metalsi&i 5 - The primary purpose of our work 
is to give an account of the application of the EMTO 
method to the calculation of ITS in bcc nonmagnetic 
metals and random alloys. Second, we aim to provide 
a consistent theoretical guide to further optimization of 
the composition of V-based alloys for structure material 
applications. 

The structure of the manuscript is as follows. In Sec- 
tion nn we describe the computation tool and all impor- 
tant numerical details. The results are presented and 
discussed in Section IIIII Here, first we assess the accu- 
racy of our calculations by considering pure Mo and V. 
Then we study the effects of the alloying elements on the 
ideal strength of Mo-Tc and V-Ti-Cr alloys. We make 
use of band filling arguments, canonical band theory, and 
structural energy differences to explain the trend of the 
computed ideal strength of V-based alloys. 



II. COMPUTATIONAL METHOD 
A. Total energy calculation 

The first-principle method used in this work is based 
on density functional theory (DFT)^ 4 - formulated within 
the Perdew-Burke-Ernzerhof (PBE) generalized gradient 
approximation for the exchange-correlation functional. 35 
The Kohn-Sham equations were solved using the EMTO 
method) 3 ^— The problem of disorder was treated within 
the CPA and the total energy is computed via the full 
charge-density technique i^2r— 

The EMTO method is an improved screened Korringa- 
Kohn-Rostoker (KKR) method^ where the full poten- 
tial is represented by overlapping muffin-tin potential 
spheres. Inside these spheres, the potential is spheri- 
cally symmetric and constant in between. By using over- 
lapping spheres, one describes more accurately the ex- 
act crystal potential compared to conventional muffin-tin 
or non-overlapping methods. Further details about the 
EMTO method and its self-consistent implementation 
can be found in previous works i^ ~ 38 i 41 ' 42 The accuracy of 
the EMTO method for the equation of state and elastic 
properties of metals and alloys has been demonstrated 
in a number of previous works j 32 i 37 i 40 i 43 ~ ilHowever , to 
our knowledge, no CPA-based alloy theory has previously 
been employed in the ab initio determination of the ideal 
strength of random solid solutions. This may be ascribed 
to the limited accuracy of the classical CPA-based elec- 
tronic structure methods for systems with reduced crys- 
tal symmetry. In this respect, the present work repre- 
sents a pioneering demonstration of the EMTO-CPA ap- 
proach in ideal strength calculations for random alloys. 



B. Ideal tensile strength calculations for bcc 
crystals 

The first step in ITS calculations is to compute the 
equilibrium lattice constant of the material in the ground 
state structure. The present alloys adopt the bcc struc- 
ture. In the second step, an uniaxial tensile strain e is 
applied along a specific crystalline direction which mim- 
ics a certain tensile stress cr. For each value of the strain, 
we relaxed the deformed structures under constraints de- 
scribed below to make sure that no internal forces remain 
in the crystal in directions perpendicular to the applied 
stress. From the above steps, we obtained energy versus 
strain curves and derived stress versus strain data. The 
first maximum on the stress-strain curve defines the ITS 
for the selected strain path. 

The stress a is given byiS 



where E is the total energy per atom and ft(e) is the vol- 
ume per atom at a given tensile strain, e is the strain of 
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c bctHct - 1 %ct/%ct > 1 <W%ct = V2 



FIG. 1. (Color online) Illustration of the lattice distortion of the bcc structure due to an applied uniaxial stress in the [001] 
direction, (a) The unit cell in the absence of strain is bcc (in the bet delineation Cbct/abct = 1). For finite strain, the lattice 
symmetry is initially distorted to bet (b) and becomes fee for Cbct/^bct = v2 (c). In panel (c), two entire fee conventional unit 
cells are drawn to ease its identification. 



the simulation cell in the direction of the applied unaxial 
force F and is defined as 



where Z|| and Iq denote the length of the cell parallel to 

F in the final state and in the initial state (without any 
force), respectively. The initial state corresponds to the 
equilibrium bcc structure with energy Eq. We define the 
uniaxial strain energy AE(l^,F) following Ref . IH as the 
total energy change upon deforming the material in the 
direction along the applied force, and relaxing with re- 
spect to the dimensions in the plane perpendicular to F, 
viz. 

AE(l«,F) = min#(ai,a 2 ,/||)-.Eo. (3) 

ai,a 2 

The minimization is done with respect to the pair of unit 
cell vectors {ai,a 2 } _L F. 

In this work, we computed the ITS in the < 001 >, 
< 111 > and < 110 > directions of bcc elemental solids 
and alloys. Three different unit cell were employed to 
describe the distorted lattices since each direction breaks 
the symmetry of the bcc structure in a different way. 
First, we chose the [001] axis to be the direction of the ap- 
plied force for < 001 > ITS calculations. The symmetry 
of the bcc lattice is reduced to the body-centered tetrag- 
onal (bet) one, see Fig. [TJ if the relaxation of the lattice 
is restricted to a change in length of the quadratic basal 
plane only. Accordingly, we may rewrite Eq. ([3]) with the 
in-plane lattice constant abet = a>i = Q>2t = |aj|, and 
ai • a 2 = 0, 

AE(c hcU [001]) = min£( abct ,c bct ) - E , (4) 

a bct 

with Cbct = in denoting the unit cell length along [001]. 
Increasing the axial ratio Cbct/flbct from 1 to y/2 trans- 
forms the bcc lattice into the fee lattice while the crystal 
remains bet during the transformation, see Fig. [T] This 
transformation corresponds to the Bain transformation 49 



and may occur under uniaxial loading environment^ Al- 
though the symmetry of the distorted lattice due to an 
unaxial load along one of the high symmetry < 001 > di- 
rections may be lower than tetragonal, e.g., orthorhom- 
bic if the relaxation is non-isotropic (di ^ a 2 )n^ the 
bet symmetry is often maintained ! 14 ' 16 Since our aim is 
to compare the present results with previously published 
ones in order to establish the accuracy of the EMTO 
method for this type of calculations, here we did not con- 
sider lower symmetry than bet for the [001] calculations. 

In an actual calculation of the stress-strain relation 
we selected a set of Cbct values equidistantly separated 
by 2.2% of the bcc equilibrium lattice parameter corre- 
sponding to different strains e. For each value of Cbct we 
computed the total energy on a dense mesh of abet values 
to find the minimum according to Eq. ([?]). The total en- 
ergies as well as the volume of each state of strain were 
fitted to polynomial fit functions and differentiated to 
obtain the stress as a function of strain (Eq. |T])). The 
maximum strain, e m , is then the ITS at the maximum 
stress, cr m . The choice of the order of the polynomial fit 
functions was crosschecked to avoid any substantial in- 
fluence on the values of er m and of e m . For consistency, 
we employed the same fit function for all types of dis- 
tortion, i.e., along the < 001 >, < 111 >, and < 110 > 
directions. 

Straining the bcc structure along the [110] axis re- 
duces the lattice symmetry to face-centered orthorhom- 
bic (fco) symmetry, assuming that relaxation preserves 
the symmetry in the (110) plane. In the absence of 
strain, the fco lattice parameters, af co , 6f co , and Cf co , ful- 
fill af co = 6f co = \/2cf co , see Fig. [2] for a detailed illustra- 
tion. Here, strain is applied to 6f co = lu while at co and 
Cf co arc relaxed keeping Cf co /af co fixed to the initial value 
of \/2 to make the computations feasible. We may hence 
write for the uniaxial strain energy of the [110] distortion, 

AE{b ico , [110]) = min£;(a fco , 6 fco ) - E Q . (5) 

atco 

If an uniaxial strain is applied along the [111] direc- 
tion parallel to the body diagonal of the bcc structure, 
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FIG. 2. (Color online) Illustration of the lattice distortion of the bcc structure due to an applied uniaxial stress in the [110] 
direction, (a) and (b) show the undistorted lattice and the distorted lattice, respectively. Both the fco unit cell used in the 
computation and the (undistorted and distorted) bcc cell are sketched. 




FIG. 3. (Color online) Illustration of the lattice distortion of the bcc structure due to an applied uniaxial stress in the [111] 
direction, (a) shows the undistorted bcc lattice with hexagonal delineation (a wedge representing one sixth of the conventional 
hex unit is drawn), the axial ratio is Chox/ihox = a/3/8. The undistorted unit cell projected along [111] ([0001]) is depicted in 
(b) to show more clearly the full symmetry of the hex lattice. A finite strain, (c), lowers the symmetry of the bcc lattice to a 
trigonal one. 



the symmetry of the lattice is reduced to trigonal symme- 
try. The distorted lattice can be equivalently described 
by a hexagonal (hex) lattice or a rhombohedral lattice. 
Here, we chose the hex lattice, as sketched in Fig. [3] 
The [0001] axis and the [1010] axis of the hex lattice are 
oriented parallel to the [111] axis and the [110] axis of 
the bcc lattice, respectively. In the absence of strain, 
the lattice parameters of the hex unit cell and the bcc 
cell are related by Ohex = V^bcc and Chcx = \/3/4abcc, 
where ahox denotes the in-plane lattice parameter of the 
hexagonal basal plane (the distance between two atoms of 
the smallest equal-sided triangle). The out-of-plane lat- 
tice parameter, Ch cx , is oriented parallel to the threefold 
stacking axis. Bcc (111) planes are stacked along [0001] 
(ABCABC stacking) and two successive planes are dis- 
placed by 1 /3 of Chex (cf . Fig. |3]) . For each applied strain, 
interatomic distances in the (0001) planes are relaxed 
while the in-plane symmetry is preserved. The uniaxial 
strain energy relaxed with respect to ai 1C x = <xi = a<i 
becomes, 

A£?(chex) [HI]) = nun.E(ai ) - E . (6) 

When the hexagonal axial ratio Chex/ohex is identical to 
\/3/2 the strained hexagonal lattice coincides with the 



simple cubic (sc) lattice. Similarly, when Chex/^hex =V3 
we recover the fee lattice. 



III. RESULTS AND DISCUSSION 
A. Equilibrium volume 

The theoretical equilibrium lattice parameter of bcc 
Mo is 3.165A, which agrees well with the experimen- 
tal data 51 3.141A. When 10 at.% Tc is added into the 
Mo matrix, the lattice constants decreases from 3.165A 
to 3.156A. This trend can be understood by the smaller 
atomic radius of hexagonal close packed Tc (1.36 A) as 
compared to that of bcc Mo (1.40 k)M 

Theoretical equilibrium lattice parameters and elastic 
properties of V and disordered V-Cr-Ti alloys were in- 
vestigated with EMTO and the CPA in our previous pa- 
per, see Ref. HH. Since our employed lattice parameters 
are identical to these earlier results, we refer the reader 
to this reference for a detailed account of the numerical 
data. Here we only restate the most important findings 
necessary to follow the ongoing discussion. 

For pure vanadium, the calculated equilibrium lattice 
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parameter is 2.998A, in good agreement with the corre- 
sponding experimental valued 3.03 A at low tempera- 
ture. For V alloyed with Ti and Cr, we found that the 
theoretical lattice constants increase (decrease) with in- 
creasing Ti (Cr) concentrations, for instance, the lattice 
constant of the alloy with 10% Cr is 2.981 A, but for the 
alloy with 10% Ti, the lattice constant is 3.02lA. Our 
theoretical results confirmed the experimental trend^ 

Our results showed that the EMTO method with the 
CPA can be used to describe the elastic properties of V- 
based alloys and gave the accuracy and available results 
in the small deformation region. Since the ideal strength 
describes the mechanical properties of the material be- 
yond the linear stress-strain response, the ideal strength 
is very important during the design process of alloys. 
Thus, based on the equilibrium lattices constants and en- 
ergies that we obtained in our former calculations^ we 
calculated the ideal tensile strength of pure V, and in- 
vestigated the ideal tensile strength of bcc nonmagnetic 
Vi-z-yCr^Tij, random alloys as a function of Cr and Ti 
(0 < x + y < 0.1) concentrations using EMTO with the 
CPA. 



B. Ideal strength of pure vanadium and 
molybdenum 

To assess reliability of our computational approach, 
we first performed the simulation of a tensile test in 
V for uniaxial loadings along the [001], [110] and [111] 
directions. The corresponding stresses as a function of 
strain are displayed in Fig. @] Along the [001] direction, 
stress increases with increasing strain up to a maximum 
of er m = 16.1 GPa at a strain of e m = 16.0%. Clearly 
different from the [001] direction are the much higher 
ideal stresses in the [111] direction and the [110] direc- 
tion where the stress reaches a maximum of 26.7 GPa at 
e m = 35.2% and 37.6 GPa at e m = 38.3%, respectively. 

In the inset of Fig. Eta), we also show stress as a func- 
tion of strain in the small deformation region (e < 0.5%). 
These strain-stress relations are linear and follow Hooke's 
law, a — E ■ e, where E is Young's modulus, which de- 
pends on the direction of the applied forced For further 
use, the explicit expression of E in terms the cubic elas- 
tic constants (Cn, C12, and C44) for an uniaxial stress in 
< 001 > directions (i?<ooi>) and in < 111 > directions 
(£<m>) are given by 



E 



<ooi> 



6C" 



B 



-E<m> — 3C44 • 



1 



1 



C44 

3B 



(7) 
(8) 



where C = (Cn - C 12 )/2 and B = (Cn + 2C 12 )/3 
are the tetragonal shear elastic constant and the bulk 
modulus, respectively. According to Fig. U{a) (inset), 
the [001] direction exhibits the largest Young's modu- 
lus, which means that the [001] direction is the strongest 



direction among the three selected ones in the small de- 
formation region, while the Young's moduli in the [110] 
and [111] directions are smaller and degenerate. Indeed, 
using the theoretical elastic parameters of bcc V— we ob- 
tain £ <O oi>=205.3 GPa and £; <m> =101.4 GPa. These 
values are in good agreement with those derived directly 
from Fig. @}a) (approximately 200 GPa and 110 GPa for 
£'<ooi> and -E<m>, respectively). This behavior at 
small strains is distinct from the one for the larger, non- 
linear deformation region, where [110] ultimately replaces 
[001] as the strongest direction. 




15 20 25 30 
Strain (%) 



FIG. 4. (Color online) The ideal strength (a) and energy (b) 
of pure V in [001], [111] and [110] directions as a function 
of strain. The insets show the trends of ideal strength and 
energy in the smaller strain region. 



The ideal tensile strength a m corresponding to the 
strain e m from our and other calculations in [001], [111], 
and [110] directions for vanadium are listed in Table HI 
First, it can be seen that for the [001] direction, projcctor- 
augmented-wave (PAW) work o 55 ' 56 i 58 reported similar 
ideal strength values, such as Nagasako et al& gave a 
value of cr m =17.8 GPa at e m =17.0% and Liu et al. re- 
ported a value of er m =19.1 GPa at e m =18.0%. Our 
EMTO results for stress and strain are slightly smaller 
than the PAW values however still in good agreement. Li 
et al~2- tabulated a value of er m =16.8 GPa at e m =16.0% 
using pseudo potentials (PP). Our results are in agree- 
ment with their values. Krenn et al. obtained an ITS 
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TABLE I. The ideal tensile strength cr m (GPa) and the corresponding strain e m (%) from our and published computations in 
[001], [111], and [110] directions for a pure vanadium crystal and for the [001] direction of molybdenum. 



element 


method 


direction/[001] 


direction/ [111] 


direction/ [110] 


<x m (GPa) 


£»(%) 


cr m (GFa) e m [7a) 


_ 1 rin n \ _ (ty\ 




hjM 1 u 


16.1 


16.0 


26.7 35.2 


37.6 38.3 




PAW Ref. 55 


19.1 


18.0 


31.0 36.0 


32.8 42.0 


V 


PAW Ref. .6 
PAW Ref. 56 
PP Ref. rg 
PP Ref. 17 


17.8 
18.9 
16.7 
12.2 


17.0 
18.2 
16.0 
8.3 








EMTO 


26.7 


12.0 






Mo 


FP-LAPW Ref. 57 
PP Ref. 15 


25.4 
28.8 


10.0 
13.0 







value of 12.2 GPa at an attainable strain of 8.3%, which 
is considerably lower compared with all other theoretical 
results. The difference may be due to their estimation of 
the ideal stress assuming a sinusoidal stress-strain rela- 
tion on the basis of Frenkel's^ and Orowan'sSS, works. 

Comparing the ITSs in [001], [111] and [110] direc- 
tions among each other, we identify the [001] direc- 
tion as the weakest one. We note that earlier ab ini- 
tio approache a 2 i 8 i 15 ' 16 i 18 ' 61 ~ — consistently identified the 
< 001 > directions as the weak directions in tension for 
elements crystallizing in the bcc structure. The ITSs of V 
along the [111] and [110] directions are clearly larger than 
the one along [001], the ideal strength along the [110] 
direction exhibiting the largest value. Cerny and Pok- 
lud a 18 ^ 63 have investigated various bcc transition metals 
and also predicted the largest stresses for the < 110 > 
directions. 

Our computed maximum stress for the [110] direction 
is larger than the only other available theoretical data, 
Ref. l55l . which may be attributed to our constraint re- 
laxation (fixed lattice parameter ratio Cf rn / af rr > = \/2). 
If this constraint is released as done in Ref. [55|, the to- 
tal energy for each lattice distortion may be lowered and 
hence the uniaxial strain-energy curve may be shallower 
compared to a constraint calculation which eventually 
lowers the ideal strength for the [110] direction. 

The ideal stress of V in the < 1 1 1 > directions is signif- 
icantly larger than the stress in the < 001 > directions. 
This may easily be understood on the basis of structural 
energy difierence o 64 i 65 (SEDs) as outlined in the follow- 
ing. The stability of a crystal structure relative to an- 
other one for one and the same element /alloy may be 
expressed in terms of their mutual energy difference, or 
often referred to as SED. We at first note that the total 
energy versus strain curve of the strain in [111] direction 
raises stronger than the other ones in the range given (the 
corresponding [001] curve even flattens), see Fig. @] The 
distorted bet lattice of the [001] deformation and the dis- 
torted hex lattice of the [111] deformation coincide with 
the fee lattice (at Cb c t/at>ct = \/2) and with the sc lattice 



(at Chox/ahcx = \/3/2), respectively. All cubic Bravais 
lattices are symmetry dictated maxima of the respective 
energy versus strain curves,— ~— i.e., the uniaxial strain 
energy must level off to the fcc-bcc structural energy dif- 
ference and to the sc-bec structural energy difference for 
an elongation along [001] and an elongation along [111], 
respectively, albeit at strains larger than the maximum 
strain (e m ). This implies a limitation on the maximum 
stress since it restricts the uniaxial strain energy (AE) 
to the structural energy differences within the strain in- 
terval (Ae) to accomplish the transformation from bcc to 
either fee or sc, i.e., the ratio AE/ Ae is bounded. We ar- 
gue that the maximum stress can be estimated from the 
previous bound on the basis of SEDs which then allows 
to understand the anisotropy of the attainable maximum 
stress with respect to uniaxial forces applied in certain 
crystallographic directions. 

For a more quantitative analysis, we below assume a 
state of strain with constant volume fixed to the theo- 
retical bcc equilibrium volume per atom, fibec- This is 
to keep the presented argument both simple and general, 
and is neither a severe approximation since the relaxed 
volumes of all three cubic structures of V, i.e., the vol- 
umes of the initial and final states of the transformation, 
differ from another by less than 6%. We then may ap- 
proximate the "true" <7 m by rewriting Eq. (fl). 



SED 



1 AE 
fibec Ae 



(9) 



where AE is either the fcc-bcc or the sc-bec structural 
energy difference, and Ae is the strain at constant vol- 
ume necessary to transform the bcc lattice into either 
the fee lattice or into the sc lattice. Readily one may 
find, Ae [001] = (1 - O.&J^/O.b 1 / 3 ~ 0.260 and Ae [m] = 
(1 - 0.25 1 / 3 )/0.25 1 / 3 « 0.587 for the [001] distortion 
and the [111] distortion, respectively. We computed for 
V AE = E {cc - E hcc = 20.48mRy/atom and AE = 

E sc - ^bec = 82.76mRy/atom. Using fl bcc = 13.4695 A 3 
and Eq. ©, the above SEDs yield c^ed _ 12 g GPa and 
ct sed _ 22.9 GPa for the [001] distortion and the [111] 
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distortion, respectively. These simple estimates are to 
be compared to our ab initio results (Table Q1 , where we 
obtained 16.1 GPa and 26.7 GPa for the respective ideal 
stresses. Thus, the anisotropy of the ideal stress may 
indeed^ be understood on the basis of SEDs, if the uni- 
axial strain energy exhibits nearby symmetry dictated 
maxima. 

The just presented argument, however, can not be used 
to explain the relative magnitude of cr m for uniaxial load- 
ing in the < 110 > directions of the bcc lattice, because 
the distorted orthorhombic lattice (Fig. [2|) does not co- 
incide with a higher symmetric one for any value of the 
strain e > 0, i.e., the uniaxial strain energy does not level 
off as if there is a nearby symmetry dictated extremum. 
It was pointed out^ that such extrema can be created if 
the transformation path includes an identical structure 
related to the original one by shear, rotation or alternate 
structures which are connected to the original one by a 
structural phase transformation. 

To further assess the accuracy of our computational 
approach, we calculated the ITS of bcc molybdenum. 
Since the [001] direction was already identified to be the 
weakest direction of bcc Mo, we here concentrated on 
this direction only. The ideal strength of solid Mo ob- 
tained in this work and previously published results are 
shown in Table |H One can see that our present results are 
very close to the other theoretical evaluations. In detail, 
the EMTO ideal strength reads 26.7 GPa at e m =12.0%. 
Joshi et al£l obtained for the stress and the correspond- 
ing strain 25.4 GPa and 10%, respectively, using the full 
potential linearized augmented plane wave method (FP- 
LAPW). 

Based on the above findings, we conclude that our the- 
oretical tool is able to describe the ideal strength of the 
V and Mo systems with sufficiently high accuracy and 
thus can be further employed to study V- and Mo- based 
alloys. 



C. Poisson ratio 

In the following structural analysis, we investigate the 
volume-strain dependence for applied stress in the [001] 
direction and relate it to the tetragonal Poisson ratio. 
We define the tetragonal, i.e., strain dependent, Poisson 
ratio, v(e) : as the negative of the derivative of the length 
normal to the applied force with respect to the length of 
the cell in the direction of the applied force, viz. 



dof dabct(e) 
dcbct(e) 



(10) 



Figure [5] depicts the strain dependence of the volume and 
the tetragonal Poisson ratio for both V and Mo. When 
e goes to zero, the present Poisson ratio should reduce 
to the cubic Poisson ratio along the [001] direction. For 
the present metals, the experimental cubic Poisson ra- 
tios are 0.258^. and 0.338^ for Mo and V, respectively. 
The cubic Poisson ratios according to Fig. [S] are 0.20 and 



0.32 for Mo and V, respectively, which are close to the 
experimental data. 

From Eq. (fTU]) . one may readily obtain the following 
relation for the stationary state of the volume change 
with strain, df2(e) / 'de|ej = 0, 



flbct( e i) 

2cbct(ei) 



^bcc(l + e<) 3 



(11) 



where v.i is the tetragonal Poisson ratio at which the vol- 
ume change is zero. In response to the applied strain, the 
volume increases up to a critical strain (q) of approxi- 
mately 15% and 20% for Mo and V, respectively. The 
volume change at the critical strain is zero, and becomes 
negative for strains larger than the critical strain. The 
critical strains of Mo and V are larger than their ideal 
strains. From Fig. [SJ we identify the critical Poisson ra- 
tios j/; (Mo) = 0.41 and z^(V) = 0.39. These critical 
Poisson ratios occur at axial ratios Cb c t/obct(Mo)=1.22 
and c bct /a bct (V)=1.28 (cf. Eq. dJ). 
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FIG. 5. (Color online) The atomic volume and the tetragonal 
Poisson ratio of V and Mo as a function of strain in the [001] 
direction. 



D. Ideal strength of molybdenum-technetium alloy 

Next, we investigated the alloying effect of technetium 
on the ideal strength of bcc molybdenum-technetium al- 
loys in the concentration range up to 10 at.% Tc along 
the [001] direction. The stress-strain dependence for Mo- 
Tc alloys is displayed in Fig. [5] and the corresponding 
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of 9.4 GPa is close to the 10 GPa obtained by Li et al. 19 
for the same alloy composition. More importantly, the 
alloying effect obtained by Li et al. is —6.7 GPa/30%Ti, 
which is in perfect agreement with the one obtained by 
us. 



6 9 
strain (%) 

FIG. 6. (Color online) The ideal tensile strengths of Mo single 
crystal and the Mo-lOTc alloy as a function of strain in the 
[001] direction. 



data are listed in Table HU From Fig. [6j we identify 
cr m = 24.9 GPa at e=11.6% for the Mo-lOTc alloy. Com- 
pared to pure Mo, Tc decreases the ideal strength. Bcc 
Mo has a d-band occupation of approximately 4.6 elec- 
trons as obtained from our EMTO calculation. Accord- 
ing to band filling theory and DFT calculations , 64 ' 69 the 
fcc-bcc SED of Tc is smaller than the one of Mo which 
means a destabilization of the bcc structure with respect 
to the fee structure. Thus, based on Eq. the ideal 
strength of Mo-Tc should be smaller than that of Mo, in 
line with the results from Fig. [SI 



E. Ideal strength of V-based alloys 

In the following we turn to the V-based alloys and in- 
vestigated the effect of the alloying elements Cr and Ti 
on the ITS of bcc V. The total concentration of both so- 
lutes was varied in the range from to 10%. Figure [7] 
shows the composition dependence of the ideal strength 
of ternary bcc V-Cr-Ti random alloys along the [001], 
[111], and [110] directions and the corresponding numer- 
ical data for selected compositions are listed in Table [TTJ 
From Fig. [3 it can be clearly seen that the ideal strengths 
of the vanadium binary alloy for all directions increase 
with increasing Cr content and decrease with Ti addi- 
tion into pure V. For instance along the [001] direction, 
an addition of 5% Ti to V reduces the ITS by 1.2 GPa, 
i.e., from 16.1 GPa to 14.9 GPa. If however 5% Cr are 
added to the V matrix, the ITS increases by the same 
amount to 17.3 GPa. 

Li et al~2- calculated the ideal strength of Ti-V alloys 
for concentrations of Ti > 30 at.% based on a PP method 
and the VCA. Their results confirm our trend, namely 
that the more Ti is present in the alloy the more the 
ideal strength is decreased. To further assess our results, 
we also calculated the ITS of the V-30Ti alloy. Our value 




FIG. 7. (Color online) The ideal strength of V-based alloys 
in [001], [111] and [110] directions as a function of Cr and Ti. 
All stress values are in GPa. 



The addition of chromium or titanium into the vana- 
dium matrix leads to opposite consequences for the ideal 
strength of binary vanadium alloys. It is apparent from 
our data, that the effect of both elements on cr m is how- 
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TABLE II. The ideal tensile strength <r m (GPa), corresponding strain e m from our calculations in [001], [111] and [110] directions 
for vanadium-based aollys and [001] for Molybdenum-10%Technetium alloy. 



composition 


direction/ [UU1J 




direction/ [111] 




direction/[llUJ 




<7 m (GPa) 




a m (bra) 


£m(7o) 


<r m (GPa) 




V-lOCr 


18.4 


16.2 


27.7 


34.3 


38.7 


37.6 


V-5Cr 


17.3 


16.1 


27.3 


34.4 


38.2 


38.0 


V 


16.1 


16.0 


26.7 


35.2 


37.6 


38.3 


V-5Cr-5Ti 


16.0 


16.0 


26.5 


35.3 


37.3 


38.2 


V-5Ti 


14.9 


15.3 


26.1 


35.2 


36.6 


38.5 


V-lOTi 


13.7 


15.2 


25.4 


35.3 


35.9 


38.3 


Mo-lOTc 


24.9 


11.6 











ever of similar magnitude. As a result of both the op- 
posite alloying effect and their similar order, er m is al- 
most unchanged if equal amounts of Cr and Ti are al- 
loyed to V. Hence, equi-composition V-Ti-Cr alloys pos- 
sess nearly the same ideal strength as pure vanadium 
(equi-composition alloys are situated along the diagonal 
in Fig. [7]). For example, the ideal strength of V-5Cr-5Ti 
is 16.0 GPa, which is very close to the 16.1 GPa for the 
pure V solid. 

Furthermore from our data in Table UH we infer that 
the [001] direction exhibits the lowest maximum strain 
value among all three directions. The ratio of the cal- 
culated maximum strain along the [111] direction to the 
calculated maximum strain along the [001] direction for 
all alloys varies from 2.11 to 2.32. That seems plausible 
since the strain at constant volume necessary to trans- 
form the bcc lattice into the sc lattice (Aejm]) is much 
larger than the strain necessary to transform the bcc lat- 
tice into the fee lattice (Aer oii), cf. Sec. IIIIB1 in fact 
their ratio is 2.260. This agreement may however be for- 
tunate. It is obvious that the "true" e m is much smaller 
than Ae, since the uniaxial strain energy levels off around 
the stable bcc state and around the fee and sc extrema. If 
this leveling off is presumably symmetric on the minimum 
and the maximum in energy then e m is roughly Ae/2 for 
the respective directions, which however preserves their 
mutual ratio (2.260). We further observe from our data 
in Table HIl that the strain maximum for the [001] dis- 
tortion slightly increases (decreases) with increasing Cr 
(Ti) amount in the binary alloy. The equi-composition 
alloy V-5Cr-5Ti has the same data as pure V. 

A direct comparison of the theoretical ITS with the en- 
gineering strength is not possible, since the latter is deter- 
mined by the dislocation network whereas the prior is a 
simple measure of the strength of the chemical bonds for 
an ideal crystal. The present calculations indicate that 
the ITS of V decreases with Ti and increases with Cr. 
The ultimate tensile strengths of V-a;Cr-(4-5)Ti (0 < x < 
15) polycrystalline alloys have been measured in the tem- 
perature range from room temperature to about 1000K 
by Chuang et al2& For all measurements, the ultimate 
tensile strength increased monotonically with increased 
Cr content, which agrees well with the present trend ob- 



tained for ITS. On the other hand, the yield strength of 
Ti-alloyed V increases with the amount of Tij 7 -^ which is 
opposite to the trend obtained for the present ITS. 

F. Electronic structure and ideal strength 

The ideal strength of materials is ultimately due to 
the atomic bonding strength. It is well established 7 -^— 
that the electronic structure and the bonding in transi- 
tion metals is governed by a narrow valence <i-band that 
hybridizes with a broader valence nearly-free-electron sp- 
band. SEDs for the transition metal series (disregarding 
the effect of magnetism in the present discussion) that 
give rise to the experimentally observed sequence of sta- 
ble crystal structures are well understood in terms of 
band filling within this band picturei 64 ! 65 ' 77 It is clear 
that mainly the gradual filling of the d-band not only 
dictates crystal structure sequences in all three transition 
metal series, but also trends of the equilibrium volume 
dependence , 76 ' 78 the cohesive energy j 7 - 7 - - — the bulk mod- 
ulus,— and elastic constants . 81 i 82 Here we argue that the 
observed trend of er m for V-Cr-Ti alloys may as well be 
understood on the basis of band filling arguments. 

Figure [8] shows the density of states (DOS) of pure 
V, vanadium alloyed with 10 % Cr and with 10 % Ti, as 
well as the DOS of the equi-composition V-5Cr-5Ti alloy 
as obtained from our CPA calculations. We notice that 
the curves of V-lOCr and of V-lOTi are almost rigidly 
shifted with respect to the one of V, while the shape of 
their DOSs is merely effected by alloying vanadium with 
adjacent elements in the periodic table. This rigid band 
shift behavior signals an increase of the d-occupation (in- 
crease of the number of valence electrons) and a decrease 
of the number of d-occupation (decrease of the number 
of valence electrons) for Cr addition and Ti addition, re- 
spectively. Thus, alloying V with Cr and Ti changes 
the band filling. The DOS of V-5Cr-5Ti is virtually not 
shifted with respect to the one of V which indicates a 
zero net gain in the number of d-electrons in the matrix. 

Since the band filling is the most important parameter 
determining structural stability in transition metals, it 
is straightforward to correlate the d-band filling due to 
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alloying to the ideal stress via SEDs in a manner simi- 
lar to our approach in Sec. IIIIBI For the refractory ele- 
ment V with a d-band occupation of approximately four 
electrons, we expect an increase of the fcc-bcc SED if 
more electrons are added into the d-band since the fcc- 
bcc SED of Cr is larger than the one of V^^S That 
is, the SEDs corresponding to intermediate band fillings 
are simply assumed to lie in between the values of adja- 
cent elements, corroborated by results of canonical band 
theory . 64 ' 83 In other words, adding electrons starting at 
the band filling corresponding to the one of V stabilizes 
the bcc structure with respect to the fee structure. The 
opposite holds for Ti alloying: the fcc-bcc SED of Ti is 
smaller than the one of V which means a destabilization 
of the bcc structure with respect to the fee structure if 
electrons are removed from the d-band. Figure [H] displays 
the uniaxial strain energy as a function of strain for V 
based alloys in [001] direction. It is apparent that with 
respect to the curve of V, the strain curve rises more 
rapidly with increasing Cr concentration and falls with 
increasing Ti concentration. Furthermore, the energy- 
strain curve of the equi-composition alloy (V-5Cr-5Ti) 
is similar to that of V. Based on the above trends and 
according to the behavior of SEDs with band filling, we 
expect that the ideal strength in the [001] direction in- 
creases (decreases) with increasing Cr (Ti) concentration. 
To the best of our knowledge, sc-bec SEDs for 3d transi- 
tion metals have not been systematically investigated so 
far apart from the FM elements Fe, Co, and Ni,— and 
Cui^ From the trend of the sc-bec SEDs we may how- 
ever deduce the behavior of cr m along the [111] direction 
as a function of band filling. 

To confirm our expectation, we computed SEDs of 
vanadium based alloys and present the results in Ta- 
ble Hill Following the same reasoning as in Sec. 1111 Bl we 
also computed cr^ ED according to Eq. ([9]) employing the 
theoretical equilibrium volumes of these bcc alloys that 
we reported prior.— Note that the influence of the volume 
change due to alloying (which enters Eq. ©) is about 2 % 
and is not primarily responsible for the observed trend 
of er m since the volume effect is much smaller than the 
effect of AE. 

TABLE III. Structural energy differences (in mRy) of the fee 
and the sc structure with respect to the stable bcc phase for 
bcc Vi~x-yCr x Ti y random alloys. Based on these energy 
differences and volumes from Ref. [33, the ideal strength o"m ED 
(in GPa) is approximated according to Eq. @. 



composition 


direction/ [001] 


direction/fill] 


Efcc E\, cc 


SED 




SED 


V-lOCr 


23.00 


14.6 


84.79 


23.9 


V-5Cr 


21.74 


13.7 


83.99 


23.5 


V 


20.48 


12.8 


82.76 


22.9 


V-5Cr-5Ti 


20.28 


12.6 


82.76 


22.8 


V-5Ti 


19.06 


11.8 


81.98 


22.4 


V-lOTi 


17.63 


10.8 


80.81 


21.9 
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FIG. 8. (Color online) The DOSs of V-based alloys show that 
alloying with a single element leads mainly to a rigid band 
shift of the DOS. The inset displays the complete valence 
band of pure V. Energies are given relative to the Fermi level 
which is indicated by a vertical dotted line. 



22 




Strain (%) 



FIG. 9. (Color online) The uniaxial strain energy of V-based 
alloys for an applied strain in the [001] direction. The higher 
the 3d-band occupation is the steeper the curve progression 
is. 



Our calculated SEDs clearly substantiate an increase 
of -Ef cc — -Ebcc with Cr addition, while the opposite effect 
is observed upon alloying with Ti. The function E sc — 
£"bcc follows the same trend. On the basis of our simple 
estimate of the ideal strength (Table llllj) , we conclude 
with respect to bulk V an increase of <7 m for the V-Cr 
binary alloys and a reduction of the ideal strength in case 
of V-Ti binary alloys. Our data affirms the larger relative 
change of er m for the [001] direction compared with the 
relative change of cr m for the [111] direction. We find 
for the equi-composition alloy V-5Cr-5Ti, that its ideal 
strength almost retains the ideal strength of pure V. 

Here we would like to bring up an interesting and im- 
portant detail related to the above correlations. For the 
[001] distortions in V-Cr-Ti alloys, it is found that the 
composition dependence of the ideal strength correlates 
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well with that of the SEDs and also with the one fol- 
lowed by E <001> , which is ultimately determined by the 
composition dependence of CM However, for the [111] 
distortions, the ideal strength vs. SED correlation is not 
followed by the £ , <m>. For instance, Cr increases the 
[111] ideal strength and the sc-bcc energy difference, but 
slightly decreases -E<m> and C44. In other words, the 
present results confirm the often quoted C vs. fcc-bcc 
energy difference correlation. 81,82 However, no such cor- 
relation seems to be valid between the C44 (or £ , <m>) 
and the sc-bcc energy difference. 

IV. CONCLUSIONS 

Using the EMTO method in combination with the 
CPA, we have investigated the ideal strength of bcc V- 
based alloys in the [001], [110], and [111] directions and 
of Mo in the [001] direction. All ideal strengths con- 
sidered here correspond to isotropic Poisson contraction, 
viz. no shape change was allowed perpendicular to the 
uniaxial stress. The present computed values for both V 
and Mo are in good agreement with previous calculations 
and thus confirm that our methodology has the accuracy 
needed for such kind of calculations. 

For Mo-Tc random alloy, we have found that adding 
Tc to Mo decreases its ITS. In the case of V-based alloys, 
we have obtained that Cr addition increases and Ti ad- 
dition decreases the ideal strength of bcc V in all three 
crystallographic directions. As a consequence, the ideal 
strength of V-Cr-Ti alloys remains virtually unchanged 
if equal amounts of Cr and Ti are introduced into the V 
host. We have shown that the observed alloying effects 
on the ideal strength can be understood on the basis of 
band filling arguments and structural energy differences. 

The strength of a material in the range of validity of 
Hooke's law is characterized by its elastic constants. In 
our previous paper^ we investigated the strengthening 
effects of Cr and Ti on the elastic constants of disordered 
V-Cr-Ti alloys and found that Cr increases the tetrago- 
nal shear modulus (C) and Ti decreases it, while alloys 
with equi-concentrations have similar C values as the el- 
emental V. The effects of Cr and Ti on the C44 elastic 
constant of V are somewhat smaller (and have opposite 
signs) compare to those obtained for C. On this ground, 



combining our findings for the elastic regime and those 
obtained in this work for the ideal strength of V-Cr-Ti 
ternary alloys, we conclude that Cr strengthens the ITS 
and C of these V-based alloys and Ti has a weakening 
effect on both ITS and C. Furthermore, it is known ) 81 i 82 
that the E{ cc — E^ cc SEDs in nonmagnetic transition met- 
als scale with their C . This relation can now be extended 
to our vanadium alloys: an increase of the <i-band fill- 
ing by alloying Cr to V leads to increasing & and ideal 
stress along the [001] direction, which scales with the cor- 
responding SED. The effect of adding Ti to the V matrix 
(decrease of the d-band filling) is opposite to that of Cr 
addition on both C" and a m . We have shown that no 
similar correlation between C44, ideal strength and sc- 
bcc structural energy difference can be established for 
the [111] mode. 

The present results offer a consistent starting point 
for further theoretical modeling of the micro-mechanical 
properties of technologically important transition metal 
alloys. Based on these achievements, we conclude that 
the EMTO-CPA approach provides an efficient and ac- 
curate theoretical tool to design the mechanical strength 
of non-magnetic bcc random solid solutions and reveal 
the composition dependence of this fundamental physical 
parameter. Nevertheless, in such applications one should 
always monitor the basic muffin-tin and single-site CPA 
errors and make sure that they remain at acceptable level 
as a function of the lattice distortion and chemical com- 
position. The extension of the present investigation to 
the case of magnetic alloys (such as the stainless steels) 
and to other Bravais lattices is in progress. 
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